Smectic phases with cubic symmetry: the splay analog of the blue phase 
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We report on a construction for smectic blue phases, which have quasi-long range smectic transla- 
tional order as well as long range cubic or hexagonal order. Our proposed structures fill space with 
a combination of minimal surface patches and cylindrical tubes. We find that for the right range of 
material parameters, the favorable saddle-splay energy of these structures can stabilize them against 
uniform layered structures. 
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Liquid crystalline blue phases exhibit true three- 
dimensional, periodic orientational order. Two of these 
phases possess cubic symmetry (BPl and BP2) while 
the third (BP3) is thought to be an isotropic melt of 
double-twist cylinders j|, [2). Recently, new phases of 
matter have been identified that possess the quasi-long 
range translational order of smectics || and, at the same 
time, three-dimensional orientational order. These three 
distinct smectic blue phases have been observed near the 
isotropic transition of these compounds: BP sm l has cu- 
bic symmetry, while BP sm 2 and BP sm 3 have hexagonal 
symmetry. The precise physical properties of these mate- 
rials have been the study of intense investigation in recent 
years [g, q|. However, there is no obvious way to incor- 
porate smectic ordering into the traditional double-twist 
tube blue phase ordering put forward by Sethna et al 
for nematic blue phases. In general, since smectic order- 
ing is incompatible with cubic symmetry, it is expected 
that any blue phase structure must include smectic de- 
fects as well as orientational defects. Though a model for 
double-twist cylinders with smectic order has been pro- 
posed j§| , the simplest variant of that model is incompat- 
ible with experimental details ||. In this letter, we pro- 
pose a new scheme for constructing a smectic blue phase 
that fills space with continuous concentric layers with 
cubic symmetry. For this construction there are elastic 
energy costs arising from non-uniform layer spacing and 
layer bending as well as condensation energy costs aris- 
ing from melted regions. However, we show that when 
the saddle-splay constant, K 2 4, is negative enough, these 
energy costs can be compensated by the gain in Gaussian 
curvature energy in the surfaces. 

The key ingredient in our construction is the obser- 
vation that saddle-splay and the Gaussian curvature are 
identical ^ for layered systems with uniform spacing. 
The saddle splay energy of a director field N is 



F ss = K 2 4 J d 3 x V • [(N ■ V) N - N ( V • N)] 



(1) 



where -K24 is a Frank constant. This term is normally ne- 
glected since it is a total derivative and thus, for most ge- 
ometries, becomes an unimportant boundary term. How- 
ever, defects introduce boundaries, and it is precisely the 
saddle-splay energy that stabilizes the standard nematic 
blue phase. Since that phase is riddled with defects, the 



saddle-splay boundary energy grows linearly with vol- 
ume. However, if space can be filled with surfaces with 
normal N then F$s has a different interpretation - if we 
replace the z integration in ([!]) with an integration over a 
Lagrangian coordinate n, which labels the surfaces, then 



F S s = -2if24 j dnj dxdy -£-K n (x,y) 



(2) 



2X 24 / dn / dxdy^g n (x, y) a n (x, y)K n (x, y) 



where a n (x,y) is spacing between layer n and n + 1 at 
(x, y), K n is the Gaussian curvature of the n th surface, g n 
is the determinant of the two-dimensional surface metric, 
and yfg~^a n = dz/dn follows from conservation of volume. 
Since the Gauss-Bonnet theorem || implies that for a 
surface of genus g the integrated Gaussian curvature is 
4tt(1 — g), we expect that higher-genus surfaces are fa- 
vored by the saddle-splay term. Note that since a n (x, y) 
is not necessarily constant, the integral (||) is not simply 
topological. Nonetheless this identification will aid us 
in our choice of smectic structures. We note that in our 
model, the saddle-splay is a measure of the layer normals, 
not the nematic director. When the nematic director fol- 
lows the layer normal these are, of course, equivalent. 
However, in type-II smectics it is possible for the saddle- 
splay of the director field to differ in its precise numerical 
value from the saddle-splay in the layers. 

We must choose our structures in the context of the 
standard bulk free energy of a smectic liquid crystal: 



F S r, 



= I d 3 J 



(v$r - 1 



2^1^ 



(3) 



where the smectic density is p cx cos (27r$/a), $(x, y, z) is 
a phase field, a is the layer spacing, B is the compression 
modulus, K\ is the bend modulus, and H = |V-N is the 
mean curvature of the layers. Though we could capture 
the same physics by replacing H with ^V 2 $ we choose 
the former for computational convenience. 

We note that there is an intrinsic frustration be- 
tween the two terms in (||). Consider a family of lay- 
ers x n (er, t) which are solutions to <I>(x n (cr, r)) = na. 
If N = V$/|V$| is the normal to x , then x„ = 
xo + naN is a solution with vanishing compression en- 
ergy: N • V$(x„) = ^y$(x + naN) = 1 and so 
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|V$| = 1. However, if k^o are the principal curvatures 
of x , then since these surfaces are uniformly spaced, 
Ki, n = Kj.o/(l + na,Ki t o) are the principal curvatures of 
x„. Thus, if the mean curvature H = h (ki + K2) van- 
ishes for some value of n, it will not vanish for any other 
value of n unless Ki^ n = 0. Conversely, if H vanishes 
everywhere and Ki^ n does not, then it is impossible to 
have equally spaced layers. Therefore it is impossible to 
have both the bending and the compression terms van- 
ish unless all the layers are flat. Thus any ground state 
must necessarily be a compromise between non-uniform 
spacing and bending. 

Though we may choose any number of initial space- 
filling surfaces we illustrate our approach with a cubic 
structure. The Schwartz P surface, or plumber's night- 
mare, is a triply-periodic minimal surface which has the 
topology we seek. Because it is minimal, H = and 
K < everywhere on this surface. Moreover, because 
of the mathematical interest in minimal surfaces, the P 
surface has a simple parametric representations which 
proves useful in computation. This surface will be the 
template within which we fill with "concentric" surfaces. 
We refer to the section of the P-surface in the unit cell as 
the P-cell and fill it in with concentric, rescaled P-cells. 
The smaller P-cells will no longer intersect the walls of 
the unit-cell. We attach cylinders to the open ends of 
the P-cells to fill in the gaps. The advantage of this 
construction is that we can easily calculate the compres- 
sion energy for the cylinders and the bending energy for 
the P-cells. Because the inside and the outside of the 
P-surface are identical, our construction self-consistently 
fills the entire unit cell and we need only focus on the 
interior of the P-cell - the remainder of the unit cell is 
filled with eight separate octants of the P-cell. Our ap- 
proach solves the frustration between compression and 
bending by building the smectic phase out of pieces that 
have no curvature (the P-cells) and pieces which can have 
uniform spacing (the cylinders) as shown in Fig. [l]. 

The energetics of this structure has four components: 
the bending energy, the core energy, the saddle-splay en- 
ergy, and the compression energy. By our construction 
there is no bending energy on the P-cells since they are 
minimal. Thus the curvature energy is nonvanishing only 
on the cylindrical sections of the unit-cell. If p is the 
radius of the largest cylinder, then geometry requires 
p a = L/4 where L is the size of the unit cell. Since, 
by construction, the length of the cylindrical region of 
radius r is 2p a — 2r and H = l/(2r), the bending energy 
is 

F B - 24^ + (4) 

where pc is the short-distance cutoff. Inside this cutoff 
region the smectic order vanishes and there is an energy 
penalty proportional to the disordered region. If the cut- 
off is comparable to the molecular scale then, in addition, 
there is a +1 nematic disclination line down the cylindri- 
cal core which contributes to the core energy as well [jfol . 




FIG. 1: Proposed P surface smectic. In (a) we show the P- 
cell and the interior structure composed of concentric P-cells 
and cylinders. In (b) we fill in the entire (cubic) unit cell with 
eight eighths of the central structure. We have cut-away one 
quarter of the volume to show interior smectic layers. 

Including the energy of the core region at the intersection 
of these lines, we have : 

64 

P core = 24(p - p c )e H p c e (5) 

where e is a line tension. The saddle-splay energy offsets 
these positive energy contributions. As in the traditional 
chiral blue phases, the saddle-splay can easily be calcu- 
lated along the defect lines. Converting the volume inte- 
gral (Q) into a surface integral we see that the cylindrical 
cores of the P-cells contribute — 2-7T per unit length, and 
each unit cell has a length 24(p — pc) of such tubes (in- 
cluding the contribution from the other P-cell in the unit 
cell). Because the integral over the P-cell in the center 
vanishes, Fss — — 487ra|if 2 4| N, where N = (p a — pc)/ a is 
the number of layers (recall that K 2 4 must be negative) . 
Like the double-twist tubes in the chiral blue phase, the 
P-cells themselves do not contribute to the saddle-splay 
though it is their topology that forces line-like topological 
defects. 

The final energetic contribution to the smectic free en- 
ergy arises from the nonuniform spacing of the P-cells. 
By our construction, the part of the smectic which con- 
sists of concentric pieces of P surface cannot have uni- 
form spacing across the layers. However, there is an en- 
ergetically preferred difference in average radius between 
consecutive layers. We can vary the layer spacing to min- 
imize the compression energy. To perform the spatial in- 
tegral we divide our proposed structure into those parts 
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which have cylindrical symmetry and those which are 
pieces of P surface. In each region we define our coordi- 
nates such that surfaces of constant phase are constants 
of one variable. In the cylindrical region we use polar 
coordinates. Since each P-cell is a minimal surface, it 
remains minimal when we uniformly rescale it. Thus for 
each concentric P-cell, we may define a new variable Q 
through 



r = (r (0,0), 



(6) 



where r, 9, and <j) are the usual spherical coordinates. In 
the central region £ runs from to 1 and so r — r Q (9, <f>) 
is the equation of the P-cell. 

If the phase field $ describes concentric layers, then 
&[(r (9, 4>), 9, 4>] must be a constant (independent of 9 
and <fi) for fixed (: therefore <&[£r (#, </>), 9, </>] can only de- 
pend on £. However, since the (£, 9, (f>) coordinate system 
is not orthogonal, |V$| can depend on all three coordi- 
nates. To calculate |V$| we note that the layer normal 
is parallel to V<I> and so r9 r $ = r ■ V$ = (r • N) |V$|. 
But since $ = $(£) = ^(r/r ), we also have rd r § = 
(r/r )<9 c $ = C<9 C $(C) and so (r • N)|V*| = C#c $ (0 and 
is thus constant on any surface of fixed £. Choosing^ an 
arbitrary reference direction (9 ,cf> ) we have (V$) = 

p(M)A(C),with 



N( 







A(C) = [V$f 



N 



(7) 
(8) 



Where p(9, <f>) is evaluated on a shell of constant $. Note 
that p(9, 4>) is scale independent and may be evaluated on 
any of the concentric P-cells. By choosing the reference 
direction to be a point on the interface between the P-cell 
region and the cylindrical region we can define £ = pj ' p a 
inside each of the six cylindrical regions, where p is the 
radius in cylindrical coordinates. This choice allows us to 
naturally continue A(£) into the cylindrical regions. The 
angular and £ integrations decouple in these coordinates: 



Fn = 



B 



C 2 (/pi-/p2A(C) + /p 3 A 2 (C)) 



+24<(W)(l-A(C)r 



(9) 



where the three numerical constants Ip\ w 72.42, Ip2 ~ 
131.48, and Ip 3 w 68.38 are moments of p(6, (f>) which 
capture the geometry of the P-cell and were calculated 
with the aid of the the Surface Evolver software pack- 
age |pH . Since the other contributions to the total free 
energy do not depend on A(C), we minimize Pc by vary- 
ing the integrand in (||) with respect to A(£) and find 
that the resulting total compression energy in the lay- 
ered P surface structure is Fc ~ 1. 54Bpl, with a very 
weak dependence on the core size pc for pc <C p - 

For our proposed smectic structure to be stable against 
the uniform flat phase, the positive energy contributions 



from Pc, Pb, and P coro must be compensated by a large 
negative saddle-splay energy. Choosing the cutoff p c to 
be the molecular scale a, we minimize the free energy 
density T = (4p )- 3 [F c + F B + F COIC + F ss ] with re- 
spect to p and find a stable minimum when 
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e/tt, 



\K M \ > Pi Nln TV 
\K U \ » Bp 2 N 

The preferred cell size is on the order of: 

37r|P 24 | + [4/^-3/2]£ 



MK 2 A 



(10) 



(11) 



It is worth commenting on the limits in (10). At the 
smectic to nematic transition, neither K\ or P 2 4 suffer 
from anomalous divergences. Thus this phase is stable 
only when | 1 is exceptionally large. Numerically min- 
imizing to give a concrete example, we find that when 
P 2 4 = -500Pi and e = 0.99(27r|P 2 4|) that A^ » 200. 
For this structure to be stable against the flat phase, Bp\ 
must be smaller than ~ 80 K\ . Since B — K\ / A 2 where A 
is the penetration length, this requires that A > 24a, i.e. 
the system must be an extreme type II smectic. Were 
we to relax the constraint that H = in the central re- 
gion, the compression and curvature could compete and 
we would expect that Pc ~ \JBK\p 2 Q ~ B\p 2 Q . In this 
case A could be significantly smaller, though it would still 
scale as N. 

We note that if this smectic phase is stable then the 
equivalent nematic phase with B = would also be sta- 
ble. This nematic texture is the splay-analog of the tra- 
ditional chiral blue phase. This splay version of the blue 
phase arises through a completely different mechanism 
than the traditional cholesteric blue phase. In the tra- 
ditional blue phase the defects were kept apart by the 
pitch of the cholesteric; here the defects are kept apart 
by the large energy cost of the line defects which connect 
them. Moreover, the traditional blue phase is stabilized 
by a large positive value of P 2 4 as opposed to the large 
negative value required here 0. Though it is likely that 
the core energy density e should be greatly reduced in 
the case of the nematic since the smectic condensation 
energy no longer contributes to the energy, it is possible 
to have a nematic phase with a periodic lattice of defects. 
However, as e shrinks, the number of layers approaches a 
limiting value of 1.5 and the lattice constant is compara- 
ble to the core size. This is a pathological limit and, if we 
were to include fluctuations, it realizes the melting of the 
nematic phase through topological defects. The remark- 
able aspect of both the nematic and the smectic phases 
presented here is that they exist without any chirality of 
the molecules. Since the stability of this phase requires 
a careful tuning of the core energy density s and P24, it 
can only occur over a small range of temperature. 

To determine the degree of ordering in our structure, 
we calculated the powder-averaged form factor of the unit 



4 



0.35 r- 
0.3 - 
■g 0.25 - 
2 0.2 - 



ns 

0.15 ■ 



S 0.1 - 

c 

0.05 - 


50 



umbilic 




cylinder 



100 



150 200 
k ( 2a (lattice units)" 1 



250 



300 



energy would be even more negative than for the P-cell. 
Calculation of the total energy of an I-WP smectic is 
more difficult than for the P smectic because the region 
of the cubic cell outside the I-WP surface is not identical 
to the region inside. Based on the results here we expect 
that the compression energy will be relatively unimpor- 
tant for sufficiently large A/a. 

We have outlined a model for smectic blue phases and 
the splay-analog of the traditional blue phase. The fa- 
vorable saddle-splay energy of our structure is sufficient 
to stabilize it against the flat smectic phase for a realis- 
tic range of material parameters which may be encoun- 
tered near the isotropic transition. The model includes 
the flexibility to fit previous data on physical smectic 



FIG. 2: Rotationally averaged form factor for the unit cell 
with N = 50. The umbilic peak corresponds to the layering 
along the diagonal of the unit cube. The cylinder peak arises 
from the concentric cylinder regions. 



cell for N = 50. Since every line passing through the cen- 
ter of the unit cell intersects precisely 200 layers, we can 
infer the degree of ordering along any direction simply 
from the location of the peak. Our results are shown 
in Fig. 0. Not surprisingly, we found a large peak in 
the directions normal to the eight points at which the 
P surface is perfectly flat (the umbilics of the P sur- 
face). The corresponding direction relative to the unit 
cell are (±1,±1,±1). We also found weaker peaks along 
the (±1,0,0), (0,±1,0) and (0,0, ±1) directions. These 
arise from the concentric cylinder regions. We note that 
the eight peaks along the cell diagonal do not correspond 
to the observed symmetry of the smectic blue phases, 
though this can be changed with a different choice for 
our unit cell. However, we have demonstrated a "proof- 
of-principle" that the construction of splay blue phases is 
possible. A better candidate for BP sm l is Schoen's I-WP 
surface, pictured in Figure 0. This surface is also triply- 
periodic, though it is flat along the (±1,0,0), (0, ±1,0) 
and (0, 0, ±1) directions and so in X-ray we would expect 
the primary peaks along these directions. Since the I-WP 
surface has genus g — 7, the integrated surface Gaussian 
curvature is —24ir and we expect that the saddle-splay 




FIG. 3: Repeat unit of Schoen's I-WP surface. This surface 
has genus g — 7 so f KdS = —24tt. 



blue phases. More importantly, our model presents an 
entirely new and promising organizational principle for 
smectic systems. 
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